arXiv:1504.01037v3 [math.AP] 6 Oct 2015 


SHARP HIGH-FREQUENCY ESTIMATES FOR THE HELMHOLTZ 
EQUATION AND APPLICATIONS TO BOUNDARY INTEGRAL 

EQUATIONS 

DEAN BASKIN, EUAN SPENCE, AND JARED WUNSCH 


Abstract. We consider three problems for the Helmholtz equation in interior and 
exterior domains in (d = 2, 3): the exterior Dirichlet-to-Neumann and Neumann- 
to-Dirichlet problems for outgoing solutions, and the interior impedance problem. 
We derive sharp estimates for solutions to these problems that, in combination, give 
bounds on the inverses of the combined-field boundary integral operators for exterior 
Helmholtz problems. 


1. Introduction 


Proving bounds on solution of the Helmholtz equation 

Am + = —/ (1) 


(where / is a given function and /c G M \ {0} is the wavenumber) has a long history. 
Nevertheless, the following problems have remained open. 

(i) Proving sharp bounds on the Dirichlet-to-Neumann (DtN) or Neumann-to-Dirichlet 
(NtD) maps for outgoing solutions of the homogeneous Helmholtz equation (i.e., 
equation ([T]) with / = 0) in exterior nontrapping domains. 

(ii) Proving sharp bounds on the solution of the interior impedance problem (IIP) for 
general domains, where this boundary value problem (BVP) consists of ([T]) posed 
in a bounded domain with the boundary condition 


du 

dn 


— - ipu = g 


( 2 ) 


where 5^ is a given function and 77 G M \ {0}. 

This paper hlls these gaps in the literature. 

The motivation for considering the exterior DtN and NtD maps for the Helmholtz 
equation is fairly clear, since these are natural objects to study in relation to scattering 
problems. The motivation for studying the HP is two-fold: 


(i) It has become a standard model problem used when designing numerical meth¬ 
ods for solving the Helmholtz equation (see Section I5T] below for further expla¬ 
nation), and to prove error estimates one needs bounds on the solution of the 
BVP. 

(ii) The integral equations used to solve the exterior Dirichlet, Neumann, and 
impedance problems can also be used to solve the HP; therefore, to prove bounds 
on the inverses of these integral operators, one needs to have bounds on the so¬ 
lution of the HP - we discuss this more in 1|6] below. 
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This paper may be regarded as a sequel to [l5] and [73] as it variously sharpens and 
generalizes estimates obtained in those works. We will refer to these papers for many 
of the basic results. Although the results proved here hold for any dimension d > 2, we 
state them only in dimensions 2 and 3, hrstly since these are the most interesting for 
applications, and secondly since this avoids re-proving background material only stated 
in the literature in these low dimensions. 


1.1. Statement of the main results. Let C d = 2, 3, be a bounded, Lipschitz 
open set with boundary T := dQ-, such that the open complement 12 + := \ 12 _ is 

connected. Let 7 + denote the trace operators from 12+ to L, let denote the normal 
derivative trace operators, and let Vr denote the surface gradient operator on L. Let 
Br := {x : |a:| < R}. 

Definition 1.1 (Nontrapping). We say that 12+ C d = 2, 3 is nontrapping ifV 
is smooth (C°°) and, given R > sup 3 ,gf^_ \x\, there exists a T{R) < 00 such that all the 
billiard trajectories (in the sense of Metrose-Sjostrand [55]^ that start in 12+ fl Br at 
time zero leave 12+ fl Br by time T{R). 

Definition 1.2 (Nontrapping polygon). If G is a polygon we say that it 

is a nontrapping polygon if (i) no three vertices are collinear, and (ii), given R > 
suPxeO- 1^1 j there exists a T{R) < 00 such that all the billiard trajectories that start in 
12+ n Br at time zero and miss the vertices leave 12+ fl Br by time T{R). (For a more 
precise statement of (ii) see mm) 

Definition 1.3 (Star-shaped). Let 12_ C M'^, d = 2,3, be a bounded, Lipschitz open 
set. 

(i) we say that 12_ is star-shaped if x ■ n{x) > 0 for every x G L for which n{x) is 
defined (where n{x) is the normal to x & T). 

(ii) we say that 12 _ is star-shaped with respect to a ball if there exists a constant c > 0 
such that X ■ n{x) > c for every x G L for which n(x) is defined. 


Theorem 1.4 (Bounds on the exterior DtN map). Let u G dd[Q(,(12+) satisfy the 
Helmholtz eguation 

Au -|- k'^u = 0 in 12+ 

for A; G M \ {0} and the Sommerfeld radiation condition 

1 


( 3 ) 


du 

— - iku = o 

or 


rid-l)/2 


( 4 ) 


as r := |x| —)■ 00 , uniformly in x := x/r. If either 12+ is nontrapping (in the sense 
Definition \l . 1\) or 12_ is a nontrapping polygon (in the sense of Definition \1.B) or 12_ 
is Lipschitz and star-shaped (in the sense of Definition \1.3\f i) ). then, given ko > 0, 

< \k\ ||7+'«lliri/2(r), (5) 

for all \k\ > ko. Furthermore, if'y+u G H^fL) then dfu G L‘^(T) and, given ko > 0, 

ll^n■*^||+2(r) ~ ll^r(7+'*^)llL2(r) + \k\ ||7+ ■*^llL2(r) (6) 


for all \k\ > ko. 
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Theorem 1.5 (Bounds on the NtD map). Let 12+ be nontrapping (in the sense 
Definition \l.l\) and let u G satisfy the Helmholtz equation ([3]) and the Som- 

merfeld radiation condition (01). Let (3 = 2/3 in the case when V has strictly positive 
curvature, and (3 = 1/3 otherwise. 

Then, given /cq > 0, 

||7+M||Hi/2(r) ^ ||^n«L-i/ 2 (r) > (7) 

for all \k\ > kf). Furthermore, if dfu G L^(T) then 7 +M G H^i/L) and, given /cq > 0, 

l|Vr(7+M)||+2(r) + \k\ ||7+M||+2(r) < \k\^~^ , ( 8 ) 

for all \k\ > kf). 

By considering the specific examples of F the unit circle (in 2-d) and the unit sphere 
(in 3-d) and using results about the asymptotics of Bessel and Hankel functions, it was 
shown in [73l Lemmas 3.10, 3.12] that the bounds (0]) and ([ 6 ]) are sharp, and that ([7]) 
and dH]) are sharp in the case of strictly positive curvature. 

We prove the DtN bound (j 6 ]) and can then get a bound on the DtN map between 
a range of Sobolev spaces by interpolation. Of this range, the bound (0]) is the most 
interesting (since it is between the natural trace spaces for solutions of the Helmholtz 
equation) and thus we state it explicitly; similarly for ([H]) and ([7]). 

Our next result concerns the IIP under the following assumption about the impedance 
parameters rj. We permit a more general assumption on rj than that specified in the 
introduction: it can be variable, and need only have nonzero real part with a linear 
rate of growth in k. 

Assumption 1.6 (A particular class of rj). ri{x) := a{x)k + ib{x) where a,b are 
real-valued C°° functions onT, b > 0 on P, and there exists an a_ > 0 such that either 

a{x) > a_ > 0 for all x eT 01 — a{x) > a_ > 0 for all x G P. 

For purposes of obtaining estimates valid down to /c = 0 (and in particular, to make 
contact with applications in the work of Epstein, Greengard, and Hagstrom |22]) we 
will also state another, stronger, set of hypotheses on rj. 

Assumption 1.7 (Another class of p). p{x) := a{x)k -\- ib{x) where a,b are real¬ 
valued C°° functions on P and there exists a_ > 0, 6 _ > 0 such that 

a{x) > a_ > 0 for all x eT and b{x) > > 0 for all x G P. 

In our discussion of the impedance problem, we use 12 to denote the domain where 
the IIP is posed (instead of 12_), since we do not need the restriction that we imposed 
on 12_ that the open complement is connected. 

Theorem 1.8 (Bounds on the solution to the interior impedance problem). 

Let Q be a bounded (7°° open set in 2- or 3-dimensions with boundary P. Given g E 
L^(P), / G L^(12), and p satisfying Assumption \1. 61 let u E be be the solution 

to the interior impedance problem 

Au -|- k'^u = —f in 12 and dnU — ip'ju = g onT. 


( 9 ) 
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Then 

for all k E ^ 
replacing \k\. 


||V«|| 


L2(D) 


\k\ 


I'*^IIl2(d) ~ 


L^{n) 


+ 


L2(r) 


(10) 


If the stronger Assumption | i. 7 holds, estimate (fTOD holds with 1 + \k\ 


The bound ffTU]) is sharp. Indeed, in [721 Lemma 4.12] it was proved that given any 
bounded Lipschitz domain, there exists an / such that the solution of the IIP with 
g = 0 and this particular / satishes |fc||| m||l 2 (q) > ||/||l 2 (q). Furthermore Lemma 1531 
shows that if hi is a ball and / = 0 then there exists a g such that the solution of the 
IIP with / = 0 and this particular g satishes |fc|||M|| 2 , 2 (f 7 ) > || 5 '||L 2 (r). 

Note that Assumption 11.61 includes the cases p = ±k, and thus the bound (fTOj) 
holds for the two most-commonly occurring impedance boundary conditions, namely 
dnU — ikpu = g and d^u + ifcyu = g. 

For our application of this result to integral equations, we state a result on the 
Dirichlet trace of the solution of the IIP. 


Corollary 1.9 (Bound on the interior impedance-to-Dirichlet map). Let kl be 

a hounded (7°° domain in 2- or 3-d with boundary P. Given f E Lf{VL), g E L^(P), and 
p satisfying Assumption \l.(A let u E be be the solution to the interior impedance 

problem ([9]). Then 


II Vr(7'*^)llL2(r) + \k\ ||7'*^llL2(r) ^ II/IIl2(q) + ||5'||L2(r) 


( 11 ) 


for all k E If the stronger Assumption |i. 7 holds, estimate fllll) holds with 1 + \k\ 
replacing \k\. 


We now state two further corollaries, which are relevant for the numerical analysis of 
hnite-element discretizations for the IIP. For simplicity, we state them for \k\ bounded 
away from zero. 


Corollary 1.10 (Bound on the inf-sup constant). Let Q be a bounded C°° domain 
in 2- or 3-d with boundary F. Given f E g E iL“^/^(F), and p satisfying 

Assumption \l .(A let u E be the solution to the interior impedance problem fl9l) . 

Then, given ko > 0, 


liVnII L2(D) + 1^1 l|■*^llL 2 (Q) ^ \k\ (^||/||(_H-i(n))/ + ||fi'||_H--i/ 2 (r) 
for all \k\ > kfj. Furthermore, 


inf sup 

1 (o) f7i co'i 


\a{u,v)\ 


> 


111 ? 


( 12 ) 


(13) 


where defined by fl62]) below, is the sesguilinear form of the variational formula¬ 
tion of the interior impedance problem, and || • is the weighted H^-norm defined 

by fl22D below. 


Corollary 1.11 (Bound on the FT^-norm). Let Q be a bounded domain in 2- or 
3-d with boundary F. Given f E F^(fl), g E and p satisfying Assumption \l .(A 
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let u G be be the solution to the interior impedance problem ({9]). Then, given 

ko > 0 , 

ll'*^llH2(r2) ~ 1^1 (^II/IIl2(o) + ll5'llHi/2(r)j (14) 

for all \k\ > ko- 


Shifting to a slightly different perspective, having proved the bound flTOll for real k it is 
natural to impose the homogeneous impedance boundary condition dnU — ip'ju = 0 and 
consider the resolvent-like operator family dehned by solving the Helmholtz equation 
with this (/^-dependent!) boundary-condition. That is, we dehne 

Ri,r,{k) : L^in) L^in) 


by 

where u is the solution to 


Ri,r,{k)f = u, 

{A + k‘^)u = f 


satisfying 

dnU — i?77M = 0. 

If T] satishes Assumption 11.61 then Rj^^{k) is well dehned when /c G M \ {0}. Meanwhile, 
the strict positivity of a implies that Rj^rjik) is well dehned and holomorphic for Imfc > 
0. We can then use a simple perturbation argument to show the existence of regions 
beneath the real axis free of poles (the equivalent of “resonances” in this compact, 
non-self-adjoint setting); if we strengthen our assumptions to strict positivity of b, this 
yields a full pole-free strip beneath the real axis, while mere nonnegativity leaves the 
possibility of a singularity at k = 0. 

The following result is stated with the stronger hypothesis and consequent pole-free 
strip. 


Theorem 1.12 (Pole-free strip beneath the real axis). The operator family 
Ri^r]{k) : L^(f2) —j- L^(H) defined as the inverse of (A -|- fc^) with boundary condi¬ 
tion dnU — i? 77 M = 0, where p satisfies Assumption is holomorphic on hnk > 0. 
Furthermore there exist an e > 0 such that Rj^riik) extends from the upper-half plane 
to a holomorphic operator family on\m.k > —e, satisfying the uniform estimate 

L2(0)^.L2(q) ^ (l + |fc|)-^ (15) 

in that region. 

1.2. Discussion of previous results related to Theorems I1.4H1.8I and I1.12L 
and high-frequency estimates for the Helmholtz equation in general. The 

main previously-existing sharp bound for one of the DtN and NtD maps is the bound 
(p|) proved when is a Lipschitz domain that is star-shaped with respect to a ball 

(in the sense of Part (ii) of Dehnition II.3p . This bound was proved by Morawetz and 
Ludwig in [59] without the smoothness requirements of the boundary explicitly stated, 
but the same techniques apply to Lipschitz domains, modulo some additional technical 
work; see [TS] Remark 3.8] and [ST] Appendix A]. The DtN bounds ([5|) and ([6|) were also 
obtained in the strictly convex case by Cardoso, Popov, and Vodev in m as well as by 
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Sjostrand [72]; see also the parametrix construction in the appendix of [75]. Non-sharp 
bounds on the DtN and NtD maps were proved in 0, Ha, and ca ; see [73], §1.2] for 
a discussion of all these results. 

Of the bounds on the IIP in the literature, the only previously-existing sharp result 
was that ffTOj) holds when O is Lipschitz and star-shaped with respect to a ball. This 
was proved in 2-d when P is piecewise smooth by Melenk [CT Proposition 8.1.4] and in 
3-d by Cummings and Feng m Theorem 1]. The technical work referred to above can 
then be used to establish the bound when P is Lipschitz (see, e.g., [29] Theorem 2.6] 
where the analogue of this bound is proved for a more general class of wavenumbers). 
By the discussion immediately after Theorem 11.81 this bound for star-shaped Lipschitz 
domains is sharp. Bounds for general Lipschitz domains with positive powers of k in 
front of both ||/||L 2 (f^) and || 5 '||L 2 (r) were obtained in [251 Theorems 3.6 and 4.7], [231 
Theorem 2.4], and [73] Theorem 1.6]; see [73] §1.2] for more discussion. 

Regarding the pole-free strip result of Theorem 11.121 the analogous result for the 
exterior impedance problem follows from the exponential decay result of [1] for the wave 
equation with damped boundary conditions (in an analogous way to how Theorem 11.121 
followed from the exponential decay in fl6T]) L Furthermore, the recent work of [6l] on 
the exterior impedance problem gives quite precise bounds on the locations of poles 
much deeper in the lower half-space than those considered here. 

A crucial ingredient in the estimates obtained in this paper is the nontrapping re¬ 
solvent estimate, which we use to solve away errors for both Dirichlet and Neumann 
exterior problems. If 12+ is nontrapping, we have for any y G C'^(r2+) 

lk (^ + < <^(1 + ( 16 ) 

(see Theorem 13.11 below for a slightly rehned formulation and generalizations). This 
result follows from a combination of two separate ingredients. By work on propagation 
of singularities for the wave equation on manifolds with boundary by Melrose [56] . 
Taylor [HI], and Melrose-Sjostrand [M], we know that solutions to the wave equation 
on nontrapping domains with compactly supported initial data become smooth for 
f S> 1. A parametrix method of Vainberg [82] or the methods of Lax-Phillips [H] can 
then be used to turn this “weak Huygens principle” into a resolvent estimate (and 
indeed to obtain a region of analyticity below the real axis for the analytic continuation 
of the cutoff resolvent). The estimate (fT6|) is known to fail, by contrast, whenever there 
are trapped orbits, by work of Ralston [65]. We mainly use the estimate flTB]) as a 
black box in our estimates below, but we do need to return to the Vainberg parametrix 
construction to prove a variant of fll6p that deals with Dirichlet data for the nontrapping 
Neumann resolvent fLemma 14.31 1 

1.3. The main ideas used to obtain Theorems H.dlfTTHl and 71.121 We now give a 
brief overview of how the main results were obtained, with more detail naturally given 

inP151 

In contrast to the proofs of the bounds on the NtD map and IIP, our proof of the DtN 
map bounds in Theorem 11.41 takes places solely in the setting of stationary scattering 
theory, i.e., we never consider the associated problem for the wave equation. We use a 
“gluing” argument, where outgoing solutions for the far-held are “glued” to solutions 
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of an “auxiliary problem” in a bounded region. This type of argument goes back at 
least to Lax and Phillips [63l §5] and was used to obtain (non-sharp) bounds on the 
DtN map in [13] and 1731 . Our contribution is to choose a different auxiliary problem 
to that considered in [13] and 1731 , with this change then yielding the sharp result. 

The main ingredient for our proof of the NtD map bounds in Theorem 11.51 is a col¬ 
lection of restriction bounds for solutions of the wave equation with Neumann bound¬ 
ary conditions due to Tataru [78]. These are used in conjunction with the Vainberg 
parametrix construction briefly discussed in 1 11.21 above. 

For the bound on the IIP in Theorem 11.81 we use the results of Bardos, Lebeau, and 
Rauch [7] on exponential decay of the energy of solutions of the wave equation with 
damped boundary conditions, with the estimate (ITOll obtained by a Fourier-transform 
argument. Once fITU]) has been established for /c G M, the pole-free strip result in 
Theorem II. 121 then follows by a standard perturbation argument. 

1.4. Application of the above results to integral equations. As mentioned above, 
the results of Theorems 11.41 11.51 and 11.81 can be applied to integral equations. Our 
main result in this direction concerns the standard integral equation used to solve the 
Helmholtz exterior Dirichlet problem. 

When u is the solution to the Helmholtz exterior Dirichlet problem, the Neumann 
trace of u, d'^u, satishes the integral equation 

= /m (17) 

on F, where the integral operator is the so-called combined-potential or combined- 
field integral operator (dehned by flTH]) below), fk,n is given in terms of the known 
Dirichlet data y+n (see fITTIB . Usually the parameter y is a real constant different from 
zero, but in fact p will also be allowed to be a function of position on F. 

We introduce the notation that denotes the exterior DtN map, as a mapping 

from for |s| < 1/2, and Pfi'fi denotes the interior impedance- 

to-Dirichlet map, as a mapping from Pl^~^P{V) —)■ for |s| <1/2 (see M 

below and [T3l Theorems 2.31 and 2.32] for details on how these maps are dehned for 
these ranges of spaces). 

The inverse of A/ ^ can be written in terms of the exterior DtN map ^^^*1 interior 
impedance to Dirichlet map Pfi’fi as follows 

i^k,Tj) ^ ~ 7 — {P^tisi — iri)Pj^fi] (18) 

this decomposition is implicit in much of the work on the combined-potential operator 
A/but (to the authors’ knowledge) was hrst written down explicitly in [TS] Theorem 
2.33]. We give another, more intuitive, proof of this result in Lemma [6.11 below. 

The operator A), ^ is usually considered as a operator from L^(F) to itself (the reasons 
for this are explained in 1|6D and the bounds on the exterior DtN map and interior 
impedance-to-Dirichlet map in Theorem 11.41 and Corollary 11.91 immediately yield the 
following bound on ||(A/,j)"i||L 2 (r)^L 2 (r). 
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Theorem 1.13. Let 12+ C d = 2,3, be a nontrapping domain and suppose that p 
satisfies Assumption M.(A Then, given ko > 0, 

ll("4'fc,^)“^|U2(p)^+2(p) < 1 (19) 

for all \k\ > ko- 

Since the proof is so short, we include it in this introduction. The spaces Hl(T) used 
below are weighted Sobolev spaces dehned in 1|2] (in particular, see equation fl2^ h 

Proof. The decomposition flT^ implies that 


\\{^'k,ri) ^||l2(P)_^l2(p) < 1 + ||Pj 


-DtN||jpi(p)_j.i2(p) 


I 

p RD||p2(p)^p^l(p) 


\v\ 


p nD||L2(r)^L2(r) ■ 


( 20 ) 


Theorem [131 implies that H-PotTvIl//! 

1 (and thus \\PuD\\L^r)^mr) ~ 1^1 
immediately give ffT^ . 


(P)^p 2 (p) ^ land CorollaryiLllimplies that |lp 2 (p)^jpi(p) 
These results, along with the assumption on rj, 

□ 


-1 


We make two immediate remarks regarding Theorem 11.131 

(1) The bound flT^ is sharp, since it was proved in [121 Theorem 4.3] that || ,^)“^||p 2 (p) > 

2 when part of T is and d = 2,3. 

(2) In this paper we focus on the direct integral equation for the exterior Dirich- 
let problem, i.e., the equation where the unknown has an immediate physical 
meaning (in this case, it is the Neuman trace dfiu) but an analogous bound 
to fIlOp holds for the inverse of the operator involved in the standard indirect 
integral equation (where the unknown of the integral equation does not have an 
immediate physical meaning); see, e.g., [131 Remark 2.24, §2.6]. 

There have been two previous upper bounds on || (24'^^,^)~^||i,2(p)^p2(p) proved in the 
literature; the bound 

ll("4fc,,,)“^||L2(p)_^p2(p) < 1 -h (21) 

when f2_ is a 2- or 3-d Lipschitz domain that is star-shaped with respect to a ball and 
?7 G M \ {0} was proved in [151 Theorem 4.3] using the Morawetz-Ludwig DtN bound 
and Melenk’s bound on the IIP, both discussed in 1 11.21 Furthermore, using non-sharp 
bounds on bound 

for ?7 G M \ {0} was proved in [731 Theorem 1.11] when either D_ is a 2- or 3-d 
nontrapping domain, or f2_ is a nontrapping polygon. 

An immediate application of the bound flTop is the following. An error analysis of 
the h-boundary element method (i.e. the Galerkin method using subspaces consisting 
of piecewise polynomials with hxed degree) applied to the equation flTTl) was conducted 
in [31]. This analysis required ||(Afc,^)~^||L 2 (r)_j.L 2 (r) < 1, and so covered the case when 
|? 7 | ~ fc and f2_ is star-shaped with respect to a ball, using the bound (l2T|) . Thanks to the 
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bound dm), however, this analysis is now valid when 17+ is nontrapping and rj satishes 
Assumption 11.61 (Note that the error analysis of the hp-boundary element method 
conducted in jl6], [52] only requires for some (3 > 0, and 

thus the bound (j 22 j) is sufficient for this analysis to be valid for nontrapping domains.) 

The bound dlH]), used in conjunction with the recent results of Galkowski-Smith and 
Galkowski [2H]) [31], on essentially the norm of almost completes the study of the 
conditioning of ^ in the high-frequency limit, i.e., the study of 

cond(Aj._,j) := ||Aj.,j||+ 2 (r)^L 2 (r)||(Aj._,j)"^||+ 2 (r)^+ 2 (r) (23) 

for k large. This study was initiated back in the 80s for the case when is a ball [IT] , 
[ 12 ] . [ 2 ], with the main question considered being how one should choose the parameter 
7] to minimize the condition number. The hrst works to consider domains other than 
balls were [12], [HI- discuss the implications of Theorem 11.131 and [3l] on the 
condition number of Aj. ^ and the choice of p in 

So far we have only discussed integral equations for the exterior Dirichlet problem. 
The case of the exterior Neumann problem is more subtle, and we refer the reader to 
1 16.2h l l^T^ where this is discussed. 

This subsection has discussed the application of the bounds of Theorems I1.4H1.8I to 
boundary integral equations for real k. 

In a different direction, the pole-free strip for the IIP in Theorem 11.121 has the fol¬ 
lowing two applications in the theory of boundary integral equations. 

(1) This result is used in [22], along with results from classical scattering theory 
and effectively the relation flTS]) . to show that Aj.^ is invertible for Imk > —6, 
for some <5 > 0, when 17+ is nontrapping and p satishes Assumption 11.71 

(2) The method of [85] for hnding Dirichlet eigenvalues of the Laplacian using 
boundary integral equations relies on the existence of a pole-free strip for both 
the interior and exterior impedance problems (see [SS] Remark 7.5]). The former 
is guaranteed by Theorem 11.121 and the latter is guaranteed by [1] . 

2. Notation and preliminaries 

Let 17- C d > 2, be a bounded, Lipschitz open set with boundary T := dfl-, 
such that the open complement 17+ := \ D_ is connected. We denote the exterior 

and interior traces by 7 +, and the exterior and interior normal-derivative traces by d^. 
The symbol y will denote a function in G))°(17+) that equals one in a neighborhood 
of 17_. Additional assumptions about the support of particular cutoffs will be stated 
explicitly. 

The symbol A denotes the (nonpositive) Laplacian and □ denotes the wave operator 

7 - A. _ 

Given a function u G G^(R'^ \ for some i?o > 0 and given A G C, we say that u 
satishes the Sommerfeld radiation condition with spectral parameter A if 

du f 1 \ 

+ - = 0 (+ri57Jj (24) 


as r := |x| —?■ 00 , uniformly in x := x/r. 
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We define the weighted norm 

II'*^IIr1(A) (25) 

(we use this notation with X either ^2+, r2_, or F; in the latter case the gradient is to 
be understood as the surface gradient Vr)- 

More generally, for s G M we let H^{X) denote the weighted Sobolev space obtained 
by interpolation and duality from the spaces of positive integer order 

H^{X) = {uE L‘^{X) : G L^X), for all |a| < m}. 

As usual (see e.g. [201 §4.4]) we may identify these spaces on manifolds with boundary 
with the quotient space 

= {« e -ffJ(R“)}/{« : fin, = 0}, 

An easy interpolation (see, e.g., [H]) shows that an equivalent norm on H^{X) for 
s > 0 is + |fcril*|l 2 , 2 , and we will use this fact freely below. 

We will also have occasion to consider domain of the self-adjoint operator (—A -f 
/c2)s/ 2, ^ith A denoting the (nonpositive) Laplacian with Neumann or Dirichlet bound¬ 
ary conditions and s > 0. We let resp. denote these respective domains; for 
negative s the spaces are dehned by duality: ISOl §5.A], we note 

that P)y(r2±) = Hl{Q) and so by interpolation we have 

//*(f]±) = P^(f]±), sG [0,1]. (26) 

The norm with no subscript attached, ll#!], will denote the norm throughout. 

The following lemma connects Sobolev regularity in space-time to weighted Sobolev 
regularity following Fourier transform. Let denote the inverse Fourier transform 
taking the time variable to frequency variable k. 

Lemma 2.1. Let I G M. be a bounded open interval. There exist Cj such that 

for every u G il"(R x X) supported in I x X. 

The proof is simply intertwining the elliptic operator {d^ + A) with the Fourier 
transform to obtain the result for a G N, followed by interpolation and duality for the 
general case. 

2.1. Preparatory results for proving Theorems 11.41 and 11.51 (the DtN and 
NtD bounds). The following interpolation result (which appears as [73l Lemma 2.3]) 
shows that the DtN bound ([5]) follows from ([6]), and the NtD bound ([^ follows from 
([8]). To state this result, we denote the DtN map in D+ by and the NtD map 

by P^tD (following the notation in [131 §2-7])- PotN dehned as a map from H^P{T) 
to H~^P{T) by standard results about the solvability of the exterior Dirichlet problem 
and the dehnition of the normal derivative, and the regularity result of Necas stated as 
Lemma [2.31 below implies that t)e extended to a map from Ff^(F) to L^(F). 

Analogous arguments hold for 
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Lemma 2.2. ( [7^ Lemma 2.3]) With PotN’ Pmo defined above, 

||Hi/2(r)-j.H-i/2(r) — l|-^-Dw||//i(r)^L2(r) 

and analogously, 

\\PNtD\\H-^/^{r)^H^/^{r) — ll-^NtD||L2(r)-i.Hi(r) • 

(Note that an analogous result holds for the interior impedance-to-Dirichlet map, 
and thus the bound in Corollary 11.91 implies a bound on this map from H~^P{T) to 
but we do not need this latter result in this paper.) 

Having reduced the problem of obtaining the DtN and NtD bounds in Theorems 11.41 
and 11.51 to the problem of obtained the bounds between the spaces H^{T) and L^(r), 
we now use the well-known fact that a Rellich-type identity can be used to bound 
the (highest order terms of the) DtN and NtD maps, modulo terms in the domain. 
The next lemma is a restatement of Necas’ result for strongly elliptic systems (see m 
§5.1.2, 5.2.1], [Sni Theorem 4.24]) applied to the specihc case of the Helmholtz equation, 
where we have kept track of the dependence of each term on k (see m Lemma 3.5] for 
details). 

Lemma 2.3 (DtN and NtD bounds in H^{T)—L'^{T) modulo terms in the do¬ 
main). With D+ and x o,s above, given f G Lpoj„p(D_|_), let u G be a solution 

to Au + k‘^u = —f. 

(i) If £ -^Hr) then dfu G Lf{T) and 

ll^n'^lL2(r) ~ ll^r(7+'*^)llL2(r) + IIx'w|IhRh+) + 11/1^2(0+) • (^7) 

(a) If dfiu G L^(r) then x+u G H^fT) and 

l|Vr(7+M)llL2(r) ^ ||^n“|lL2(r) + 1^1^ \h+u\\l2(^r) + IIxm||^i(o+) + II/IIl2(h+) • (28) 

Therefore, to prove the bounds in Theorem 11.41 it is sufficient to prove that, if 
u G i/j)jj,(D+) is the solution to the exterior Dirichlet problem for the homogeneous 
Helmholtz equation, with iL^-Dirichlet boundary data go, then 

\\xu\\Hiin+) ^ Il7+w|lr/i(r)- 

Similarly, to prove the bounds in Theorem [L5] it is sufficient to prove that, u G 
is the solution to the exterior Neumann problem for the homogeneous Helmholtz equa¬ 
tion, with L^-Neumann boundary data g^, then with fi as in Theorem 11.51 

h+uh^D ^ ll^n“llL 2 (r) and \\xu\\Hl(n+) ^ ll^n“|L 2 (r) 

(we will actually prove the stronger result that the second bound holds with a smaller 
power of k on the RHS, but this will not affect the bound on the NtD map). The 
asymmetry between what we need to prove for the Neumann problem versus what we 
need to prove for the Dirichlet problem is due to the fact that only the iL^-semi norm 
of the Dirichlet trace is controlled in (1281) . which is due to the structure of the Rellich 
identity (see, e.g., [73l Equation 3.13]). 

Finally, in our proof of the NtD estimates we will need the following lemma. It 
is perhaps easiest to state this in terms of norms of u over Qr := D+ H Br, where 
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Br := {x : |x| < -R}, but the result could be translated into norms of xu over 17+ for 
appropriate cut-off functions x- 

Lemma 2.4 (Bounding the norm via the norm and the data). Given 
f G L^Qj^p(17+), letu G be a solution of the Helmholtz equation Au-l-fc^M = —/ 

in 17+. Then, given R > sup^gQ_ |x|, 

l|V«lly,n.) < {kf ||«|| l2(Oh+i) + (^) ll/lli,2(Qg) + 117+"“ 11 1,2 (r) ||dj;-u||^2(r) 

for all fc G M. 

This result when one of 7 M and dnU is zero is proved in [73l Lemma 2.2]; a similar 
result appears in [58l Lemma 1]. 

3. Exterior Dirichlet-to-Neumann estimates 

In this section we prove Theorem 11.41 i.e., a bound on the exterior Dirichlet-to- 
Neumann map for solutions of the Helmholtz equation satisfying the Sommerfeld radi¬ 
ation condition. 

The methods used here will be completely in the setting of stationary scattering 
theory, i.e., we will never have recourse to energy estimates for solutions to the wave 
equation (which is, of course, connected via Fourier transform). The energy estimates 
that we present are more widely known in this latter setting, however—cf. Hormander 
m §24.1] as well as the more general estimates of Kreiss and Sakamoto in the context of 
general hyperbolic systems with a boundary condition satisfying the uniform Lopatinski 
condition [10], [SB], [Bl]- (In contrast, when dealing with the Neumann-to-Dirichlet 
operator below, we need to use results known only in the wave equation setting.) 

More specihcally, the method we use to prove Theorem 11.41 consists of a “gluing” 
argument, where outgoing solutions for the far-held are “glued” to solutions of an 
“auxiliary problem” in a bounded region; this type of argument goes back at least to 
Lax and Phillips [631 §5]. In our situation, estimates for the DtN map for a lower-order 
“perturbation” of the Helmholtz equation are used in conjunction with the resolvent 
estimate for the problem with homogeneous boundary conditions. This argument was 
hrst used to obtain bounds on the DtN map in [12], and later rehned in [72]. Both 
these previous works use the equation Aw — kfw = 0 as the lower-order perturbation, 
and obtain non-sharp bounds on the Helmholtz DtN map. Here we use the equation 
Aw + (fc^ -|- i|/c|)tc = 0 as the lower-order perturbation (i.e., the Helmholtz equation 
with some absorption/damping), and this change is sufficient to prove the sharp result. 

Before we begin, it is helpful to recall the following resolvent estimates for the Dirich- 
let problem (all but one of which hold for the Neumann problem as well). 

Theorem 3.1 (Resolvent estimates). Let f G L^(f7+) have compact support, and let 
u G iL[Q^(f7+) be a solution to the Helmholtz equation Au + lf^u = —/ in 17+ that satisfies 
the Sommerfeld radiation condition (]4|) (with X = k) and the boundary condition y+u = 
0. If either 

(a) 17+ is a 2- or 3-d nontrapping domain (in the sense of Definition \l.l\) or 

(b) 17- is a nontrapping polygon (in the sense of Definition \1.2\) . or 
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(c) f2_ is a 2- or 3-d Lipschitz domain that is star-shaped (in the sense of Definition 

\^\(^)) 

then, given ko > 0, 

\\xM\Hl{n+) ~ ll/lli,2(o+) (29) 

for all \k\ > kf). 

Proof. The result for Part (a) is proved in [82l Theorem 7] using the propagation of 
singularities results of [5l], [55]. (See also Vainberg’s book [83] for a broader survey of 
these methods.) The result for Part (b) was proved when is a nontrapping polygon 
in [21 Corollary 3]. The bound fl29p was proved when is a star-shaped domain in 2- 
or 3-d in [T^l Lemma 3.8]. □ 


Lemma 3.2. If w satisfies 

Ata-|--|-i|/c|)ta = 0 in (30) 

and the Sommerfeld radiation condition (jl]) with spectral parameter ^^k'^ i\k\, then, 

given ko > 0, 

lkllHi(H+) ^ 1^1 ll7+'«^llL2(r) Pn^Wmr) ■ (31) 

for all \k\ > ko- 


Proof. Given ko > 0, there exists a c > 0 such that Im > c; therefore, since 

w satishes the Sommerfeld radiation condition and the associated asymptotic expansion 
(see, e.g., m Theorem 3.6]), w decays exponentially at inhnity; hence both w and Vtc 
are both in L^(f2_|_). 

We can therefore apply Green’s identity (i.e., multiply the PDE fl30|) by w and inte¬ 
grate by parts), and obtain that 


'y+wdfw + (fc^ i 


w\ 


iVwp 


0 . 


Jr Jn+ 

Taking the imaginary part of this last expression and using the Cauchy-Schwarz in¬ 
equality yields 


1^1 lkllL2(H+) < ll7+'«^|lL2(r) ll^n w^|L2(r) • (32) 

Taking the real part yields 

< /c^ + ||7+M;|]^2(r) , (33) 

and combining (1321) and (1331) yields the result (l3T|) . □ 


Lemma 3.3 (Bound on the exterior Dirichlet problem with damping). 

Given go G i/^(P), let w he the solution of 

Aw-\-{k'^-\-i\k\)w = 0 in ■j^w = go on P, 

satisfying the Sommerfeld radiation condition (j4]) (note that the existence of a unigue 
solution to this problem follows from Remark \3.4\ below). Then 

II^IL^(n+) $ ll5Dll//i(r) ■ 


(34) 
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Remark 3.4 (Existence of outgoing solutions to the Dirichlet problem with 
damping). 

If w satisfies Aw + {kfi + i\k\)w = 0, then w satisfies the Helmholtz equation Aw + 
X^w = 0 with A = \/k‘^ + i\k\. Since ImA > 0, the existence of outgoing solutions 
(i.e. solutions satisfying the Sommerfeld radiation condition (jlj) j to the Dirichlet and 
Neumann problems for this equation follows in the same way as in the case ImA = 0. 
Indeed uniqueness is proved for ImA > 0 in [161 Theorem 3.13]. Existence in the case 
ImA = 0 is proved using integral equation results in [131 Corollary 2.28] (see also [T3l 
Theorem 2.10]j, but the proof goes through in the exactly the same way when ImA > 0. 

Proof of Lemma \3.3[ (i) Using the bound flTI]) in the Necas result fl271) (with w = u, f = 
i\k\w) we find that 

||-9+w||^2(r) < ||Vr(7+'«^)||i2(r) + \k\ ||7+M;|]^2(r) Pn^L^^r) ■ 
and so, absorbing the Neumann data term on the LHS we have 

||^nW^IL2(r) < l|Vr(7+'«^)llL2(r) + 1^1 ll7+'«^llL2(r) • 

Using this last expression in fIM]) . we obtain fl5T|) . □ 


Theorem 3.5 (Bounds on solutions of the Helmholtz Dirichlet problem). 

Given go G i/^(r), let u be the solution of 

Au + k‘^u = 0 in 7 +m = gn, (35) 

satisfying the Sommerfeld radiation condition (jlD (with X = k). If satisfies one of 
the conditions (a), (b), and (c) in Theorem \3.1\ then 

IIa'*^IIr1(o+) ~ WdnllnKr) ■ (36) 

Proof. Let w be as in Lemma [331 Let y G C“(f2+) be equal to one in a neighborhood 
of f2_, and define v hy v := u — yw;. This definition implies that v G and 

satisfies the Sommerfeld radiation condition (|4]) (with X = k), 

Av + k‘^v = h, and 7 +n = 0 , 


where 


h := i\k\xw — wAy — 2Vw ■ Vy. 

Since h has compact support, the resolvent estimate (l2^ implies that 

\\X'v\\Hi(n+) ~ II^IIr1(o+) > 

and thus 

\\xu\\Hi{n+) ~ II^IIr1(o+) • 

Using the bound (l34|) . we obtain the result fl36|) . 


□ 


Corollary 3.6. If satisfies one of the conditions (a), (b), and (c) in Theorem \3.1\ 
and u is the outgoing solution to the Dirichlet problem 0351) then 
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Proof. This follows from combining the bonnd fl36l) with Lemma 12.31 


□ 


Proof of Theorem l.f The bonnd ([6]) is proved in Corollary 13.61 above. The bonnd 


then follows by Lemma [2.21 


□ 


4. Exterior Neumann-to-Dirichlet estimates 

In this section we prove Theorem 11.5[ i.e. a bonnd on the exterior Nenmann-to- 
Dirichlet map for solntions of the Helmholtz eqnation satisfying the Sommerfeld radi¬ 
ation condition. 

This problem is snbtler than obtaining bonnds on the Dirichlet-to-Nenmann map, 
since the Nenmann bonndary condition does not satisfy the nniform Lopatinski con¬ 
dition, hence the classic estimates of Kreiss and Sakamoto do not apply to the wave 
eqnation, nor does the simple stationary argnment nsed above for the Dirichlet problem. 
Indeed, the problem becomes an intrinsically microlocal one, with the degeneracy of the 
normal derivative at the glancing set making even global energy estimates extremely 
sensitive to the bonndary geometry (which was irrelevant to energy estimates in the 
Dirichlet case). 

The main technical ingredient in onr argnment is a collection of estimates proved by 
Tataru [78] for solntions to the wave eqnation with Nenmann (or indeed many other) 
bonndary conditions, which we now recall. The following is a restatement of part of 
Theorem 9 of [78] . 


Theorem 4.1 (Tataru). Let T be smooth. Suppose v satisfies 

□n = 0 on D+ x [0,T], 

d+v = g, (38) 

n(0) = nt(0) = 0. 


Assume g G L‘^{T x [0,T]). Then 

V G i/"(D+ X [0,T]) 


and 

wher^ 


7+nGLf^(rx [0,T]), 
a = 2/3, f3 = 1/3 in general, 

a = 5/6, (3 = 2/3 ifT has strictly positive curvature. 


(39) 


Other results from [TS] that we shall use (Theorems 3,5) estimate Dirichlet data for 
solutions of the Helmholtz equation with homogeneous Neumann condition and interior 
inhomogeneity: 


^The positive curvature used here in dimensions d 
mental form, in general dimension. 


2, 3 generalizes to be positive second funda- 
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Theorem 4.2 (Tataru). Let T be smooth. Suppose v G satisfies 

□w = F on X [0,T], 

dfv = 0, (40) 

t;(0) = Ui(0) = 0. 

Assume F G x [0,T]). Then 

7+uGif“(rx [0,T]), 

where a is given by fl3^ . 

We now turn to an estimate analogous to the usual nontrapping resolvent estimate 
that will allow us to estimate the Dirichlet data of the Neumann resolvent for a non¬ 
trapping obstacle. 

Lemma 4.3. Assume that ^2+ is nontrapping. Let RN{k) denote the outgoing Neumann 


resolvent on (2+, 

acting on f G T^Nk- Then for /c 3> 1, for every s G M 



llx^7v(fc)x/llc“+i < \\f\\v%, 

N,k N,k 

(41) 

and for s G [0,1] 


(42) 


where a is given by fl5^ . 


We remark that 2a = 1 + (3. 

Proof. The first part of this estimate is essentially the standard nontrapping resolvent 
estimate, albeit considered in more general weighted spaces than L^. The second part 
by contrast requires Tataru’s boundary estimates together with an examination of the 
details of the Vainberg construction of a parametrix for the nontrapping resolvent [83[ 
Chapter X]. This parametrix is indeed one of the usual routes to obtaining the stan¬ 
dard resolvent estimate f flTT]) with s = 0) from the weak Huygens principle (eventual 
escape of singularities), and depends crucially on propagation of singularities results 
that enable us to conclude weak Huygens from nontrapping of billiard trajectories. For 
details, we refer the reader to Theorem 2 in [SB], Chapter X; see also [33] and [33] for 
the geometry and microlocal analysis aspects. 

To establish the first part of the result, we recall that Vainberg’s estimate (see also 
the “black-box” presentation of Vainberg’s method in m yields 

||Xii?jv(/c)X2|lL2^i2 < (43) 

We must extend to more general spaces in the domain and range. First, note that 
if (A -|- k‘^)u = —f and / has compact support in a fixed region and u satisfies the 
radiation condition then we of course can write, for xo compactly supported, 

IIxoAmII < k'^WxouW + WxofW < {k)\\f\\, 

hence for any y with smaller support than xo, 


\\au\\t) 2 < {k)\\f\\, 

N,k 


(44) 
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i.e., in particular 

\\xRN{k)x\\L2^v%^^< {k). (45) 

Thus we obtain by interpolating fH3|) and flT5|) 

\\xRN{k)x\\L2^v},^^ < 1. 

Now once again if (A + k‘^)u = —f then (A + k‘^){—A + k'^Yu = —(—A + k'^Yf, hence 
by compact support of / the resolvent estimate yields 

||xo(-A + *:")L|Lj^^<||(-A + A.y/||, 

SO that for y with smaller support than xo we have 


llx«ll<« < ll/llp«.. 

Interpolation now yields 

\\xRN{k)x\\x>^^^^x>‘+i ~ 1 

for all s > 0. Now duality (which exchanges k and —k) yields the estimate for s < 0 as 
well. This completes the proof of flTTl) . 

To prove flT2D we begin by using the Vainberg parametrix construction as presented 
in CTl to establish the estimate for s = 0. In the notation of that paper, we have (see 
the two displayed equations preceding (3.5)) 

RN{k)x = RKkYl + K{k))-^ 

where K{k) is a holomorphic family of operators that is shown to have small —)■ 

operator norm for fc ^ 1, so that (/ + K{k)) is invertible there. The parametrix R^k) 
is dehned by 

R^k) = R{k) - - Xc)Va{t)) (46) 

where Xc = 1 near r2_, and 


R{k) = -tRt^kiCHmit)x)- 


Here x (also called Xa in [T7j) is a cutoff equal to 1 in a neighborhood of hi-, H{t) is 
the Heaviside function, 


U{t) 


sint\/—A 


(sine propagator for the Neumann Laplacian), and C is a cutoff with 





t < \z\ 

h 

IV 

IT 


for some Tq > Tq. The term Va{t) is obtained by solving the free wave equation (i.e. 
with the obstacle removed) with forcing given by the error term — [□, QU{t)x and zero 
Cauchy data. Happily, its analysis will be of no concern here, as the factor (1 — Xc) 
ensures that the corresponding term in fH6|) vanishes on T. 

It thus suffices from (|46D to know that x+R{k) satishes the desired estimates. To see 
this, note that if / G CU{t)f lies in L°°([0,T]; for each T < oo, simply by the 
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functional calculus for the Neumann Laplacian and the identihcation of with 

Vlf. Now Theorem 14.21 implies that 

^+CHit)U{t)xf e if“(R X T). 

(Note that ( has compact support in time in a neighborhood of the obstacle, so there 
is no difference between local and global results here; note also that the factor of H{t) 
does not affect the regularity since f/(0) = 0.) We may now Fourier transform this 
estimate by Lemma 12.11 to get 

x+mf G m 

when f & L^. 

Finally, we extend to more general s in the estimate (|42|) . Fix Fermi normal co¬ 
ordinates near F with x denoting the normal variable (distance to F) and y denoting 
coordinates along F. Let V denote any smooth, compactly supported vector held on 
12+ such that near F, V is of the form ^ aj(x, Then V can be restricted to F to 

give a (indeed, any arbitrary) vector held Vp- Note that [A,l/] is then a second order 
diherential operator in the dy^^s only near F, hence we have (cf. [79l p.407]) 

lA U : ^ Ll 

Now if 

{A + e)u = -feHl{n+) 

with u outgoing and / compactly supported in some hxed set, then we compnte 

1/(A + k^)u = -Vf, 

hence 

(A + k‘^)Vu + [V, A]u = -Vf. 

Thus, applying the Neumann resolvent and restricting gives 

W7+W = 'y+Vu = -x+RN{k)Vf - x+RN{k)[V, A]u. (47) 

Now by the estimate fjT2l) for s = 0 obtained above, we have 

||7+7?v(fc)^/|lH.<||^/IL. <11/11^1. 

Moreover, (ITT|) yields u G ^ with norm estimated by H/Hui, hence 

Thus, again by the s = 0 estimate (l42|) . 

h+RN{k)[V,A]u\\jj^ < 

and putting together our estimate for the two terms on the RHS of flTTl) . we have 
obtained for any vector held Vp on F, 

\\Vpx^RM{k)f\\^. < ll/ll^i. (48) 

Also, just the fact that f & Lf and the s = 0 estimate gives 

(fc)||7+Riv(A;)/||^.<(fc)||/IL. <11/11^1. 


(49) 
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Since Vr was arbitrary, putting together fj^ and fj4^ yields, for / compactly supported 
in a fixed set. 


Interpolating with the s 


ll7+^iv(A:)/|| 




1 + a 


< 




0 estimate now yields fH2l) for the whole range s G [0,1]. □ 


Theorem 4.4. Let 17+ be nontrapping. For each x ^ there exists ko so that 

solutions u of the Helmholtz equation 


(A + k‘^)u = 0 in 17+ 
“|r = 9n 


(50) 


satisfying the Sommerfeld radiation condition (jl]) enjoy the bounds 
and 

for k > fco- Here a and (3 are again given by equation fl5U]) . 


Proof. Fix a cutoff function ip(t) compactly supported in (0,1) with J ip = 1. Suppose 
that Vk is the solution of 


= 0 , 

= (p{t)e~^'"*gN{y) = h^{t,y), 
V = 0 for 7 < 0. 


Note that ||hK||+ 2 (][jxo^) ^ llfi'N||+ 2 (r) for ah k; this estimate and all those that follow 
have implicit constants that are, crucially, uniform in k. 

Let I C M be an open interval containing supp (p. By Tataru’s estimates in Theo¬ 
rem ITT] (and the compact support of on I x 17+) we obtain 

II'^k|Ih“(/xo+) ~ II^K|lL2(/xr) ~ Il5'w||+2(r). 

We further choose fjit) a cutoff function supported in I and equal to 1 on supp (p. Then 
we also have 

||'0'r'K||H“(KxH+) ~ llfi'N||+2(r). 

Hence by Lemma [2.11 

Now since satishes the wave equation we have 

□ (V'u^) = [□,V']W e H^-\R X 17+) ni/"-^(M;L2(17+)) 

with the norm of the RHS again estimated by a multiple of ||5'Ar||. (Note also that ifjVK 
has compact support in 17+.) Hence sinc^ (fc)~“''’^L^(17+) C 17^“^(17+) we have 

(A + = e, G I7“-1(17+), (51) 

where 

l|eK|li,--i < IIs'aII- 

N,k 

^We are of course using the fact that a — 1 < 0 here. 
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Now the nontrapping estimates for the Nenmann resolvent as stated in Lemma [4.3Jt ell 
us that if RN{k) denotes the outgoing Neumann resolvent, then for fc 3> 1 we havqj 

IIx-Rn(^)[C fc]||//a ^ 

~ \\9n\\- 

Now consider 

u = - RN{k)[e^]. (52) 

By the foregoing discussion we have 

\\xu\\h^ < llf/ivll- 

On the other hand, we have 

(A + k‘^)u = 0, 

by construction. Moreover, since we used the Neumann resolvent in constructing u, 

= (p{k - K)gN- 

Hence if we set k = fc we obtain u as the (unique) solution of (l50|) satisfying the 
radiation condition, and have obtained the desired interior estimate. 

To derive the boundary estimates, we use Lemma 14.31 as well as Theorem 14.11 The 
latter implies that 

e Hi 

hence by fl52]) it suffices to consider the term Rjq{k)[et^. Returning to the dehnition 
flSTl) of we note that we can in fact write 

e. = + fl), where fl G i/“(/ x 0+). 

Thus we obtain a slightly rehned estimate on e^: 

e, e {k)Hnn+). 

Now since a G (0,1), the estimate fl42|) of Lemma [4.31 yields an estimate on 

as desired. (Recall that 2a — 1 = (3.) □ 

Corollary 4.5. With notation as above, 

\\xu\\Hi(n^) ^ l^r"“ll^vllL2, \k\ > ko. 

Proof. This follows from combining the bounds in Theorem l4.4l with the result of Lemma 

EH □ 

Corollary 4.6. With notation as above, we have 

II7+m||^^i < \k\^~l\9N\\L2, \k\ > ko. 

^We are using the identification of Neumann domains and Sobolev spaces for exponents in [0,1]. 
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Proof. By the second part of Lemma 12.31 we have 

h+uWnl ^ II^nII + k\\-f+u\\ + WxruWhi, 

hence the results follows from the estimates on the second and third terms given above 
in Theorem 14.41 and Corollary 14.51 respectively. □ 

Proof of Theorem li.51 The bound ([8]) follows from combining the bounds in Corollaries 
14.51 and 14.61 with Lemma 12.31 (note that 1 — /9>1 — ain both the general and positive 
curvature cases). The bound ([7]) then follows by Lemma lT2l □ 

5. The interior impedance problem 

5.1. Motivation. For readers unfamiliar with the numerical analysis literature on the 
Helmholtz equation, we explain in this section why the interior impedance problem is 
of interest to numerical analysts (independent from the fundamental role it plays in the 
theory of integral equations for exterior problems, which we discuss in ^1.41 and ^H])- 

The majority of research effort concerning numerical methods for Helmholtz problems 
is focused on solving scattering/exterior problems in 2- or 3-d (such as the exterior 
Dirichlet and Neumann problems considered in l|3]and 1|4]). Boundary integral equations 
(BIEs) are in many ways ideal for this task, since they reduce a d-dimensional problem 
on an unbounded domain to a (d — l)-dimensional problem on a bounded domain. 
However there is still a very large interest in domain-based (as opposed to boundary- 
based) methods such as the hnite element method, partly because these are usually 
much easier to implement than BIEs and partly because these domain-based methods 
usually generalize to the case when k is variable (as occurs, for example, in seismic- 
imaging applications). 

When solving scattering problems with domain-based methods, one must come to 
grips with unbounded nature of the domain. This is normally done by truncating the 
domain: one chooses a (large) bounded domain 12 D 12_, imposes a boundary condition 
on dVl, and then solves the BVP in 12 \ 12_. If 12 is a ball, one can choose the boundary 
condition on c212 such that the solution to the BVP in 12 \ 12_ is precisely the restriction 
of the solution to the scattering problem—one does this by using the explicit expression 
for the solution of the Helmholtz equation in the exterior of a ball, and the relevant 
boundary condition on c212 involves the so-called Dirichlet-to-Neumann operator (see, 
e.g., [3^ §3.2] for more details). Alternatively one can impose approximate boundary 
conditions (often called absorbing boundary eonditions or non-reflecting boundary con¬ 
ditions since their goal is to absorb any waves hitting (912 instead of reflecting them 
back into 12), the simplest such one being du/dn — iku = 0 on c212. This can be viewed 
this as an approximation to the radiation condition (HI) . 

Therefore, in the simplest case, truncating a Helmholtz BVP in an unbounded domain 
yields a BVP for the Helmholtz equation in the annulus-like region 12 \ 12_, with an 
impedance boundary condition on c212, and either a Dirichlet or Neumann boundary 
condition on P. Without a fc-explicit bound on the solution of this BVP, a fully k- 
explicit analysis of any numerical method is impossible, and therefore the problem of 
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finding /c-explicit bounds on the solution of this truncated problem was considered in 

pg, m]. 

Going one step further, although the geometry of the scatterer plays an important 
role in determining the behaviour of the solution, many features of numerical methods 
for the Helmholtz equation (such as whether the so-called pollution effect occurs) can be 
investigated without the presence of a scatterer at all; this then leads to considering the 
Helmholtz equation posed in a bounded domain with an impedance boundary condition, 
i.e., the HP (and the impedance boundary condition can then be viewed as a way of 
ensuring that the solution of the BVP is unique for all k). The problem of hnding 
/c-explicit bounds on the solution of the HP was therefore considered in ra. lai. El, 

Pg, and |73|. 

Midway between, in some sense, the truncated scattering problem and the HP are 
BVPs posed on bounded domains, where impedance boundary conditions (or more 
sophisticated absorbing boundary conditions) are posed on part of the boundary, and 
Dirichlet or Neumann boundary conditions are posed on the rest. The most commonly- 
studied such problem is the Helmholtz equation in a rectangle with impedance boundary 
conditions on one side and Dirichlet boundary conditions on the other three, motivated 
by the physical problem of scattering by a half plane with a rectangular indent (or 
“cavity”). Bounds on this problem were obtained in [H] and jT^], and the recent paper 
[TU] seeks to determine the optimal dependence on k via numerical experiments. 

5.2. Interior impedance estimates. We begin with a result about uniqueness of 
solutions of the HP for complex values of the spectral parameter k. 

Lemma 5.1 (Uniqueness of the IIP). Consider the IIP ([2]) with 

r}{x) = a{x)k -|- ib{x), (53) 

where a,h are real-valued C°° functions on P. 

(i) If there exists an a- > 0 such that 

a{x) > a_ > 0 for all x G P, (54) 

and b{x) > 0 on P, then the solution of the IIP is unique for all k ^ 0 with 
Im k >0. 

(ii) If there exists an a- > 0 such that (1541) holds and there also exists a 6_ > 0 
such that 

b{x) > b_ > 0 for all x G P, (55) 

then the solution of the IIP is unique for all k with Imk >0 (i.e. we now also 
have uniqueness when k = t)). 

Proof. If u is the solution of the homogeneous HP (i.e. / = 0 and S' = 0) then applying 
Green’s identity and using the impedance boundary condition we hnd that 

ik f a|7n|^ — f b\'yu\‘^ — ( |Vn|^ -I k^ f |n|^ = 0. 

Jr Jr Jn Jn 
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Therefore, taking real and imaginary parts, and writing k = kR + ikj with kji, kj G M, 
we have 

— kj f a\'yu\‘^ — f b\'yu\‘^ — f \Vu\'^ + {k\ — k'^) / |m|^ = 0, (56) 

Jr Jr Jn Jn 

and 

kn I a|7M|^ + ‘ikftki j |m|^ = 0 (57) 

Jr Jn 

Proof of (i): if /cj? 7 ^ 0 and kj > 0, then using the assumption flHTll on a in fl571) we see 
that 'yu = 0. The impedance boundary condition then implies that dnU = 0, and thus 
Green’s integral representation (see, e.g., [501 Theorem 7.5]) implies that n = 0 in G. 
If /c/j = 0 and kj > 0, then using both the assumption fl5T|) on a and the assumption 
that b is non-negative in fl56|) . we see that u = 0 in G. 

Proof of (ii): from Part (i) we only need to consider the case when k = 0. Using the 
assumption (l55il in (l56il . we see that yu = 0 on P, and then u = 0 in G follows from 
the steps above. □ 

We now prove Theorem 1 1.8 1 bv employing the estimates of Bardos-Lebeau-Rauch [7] 
for the wave equation with the damping boundary condition, i.e. 

□n = 0 on G, (58a) 

{dn + ajdt + hy)v = 0 on P (58b) 

where a, b are smooth, real-valued functions on P with a strictly positive and b nonneg¬ 
ative. 

First we give a short proof of the standard energy estimate for the wave equation, but 
now considering the boundary condition (l58l) instead of the usual Dirichlet or Neumann 
ones. 

Lemma 5.2. Let F E L^(M x G) and G E L^(M x P) be supported in t > 0 and let v 
solve 

□n = F on Q, 

{dn + ajdt + b'y)v = G on P, 

n = 0 for t < 0 , 

where a, b are smooth, real-valued functions on P with a strictly positive and b nonneg¬ 
ative. Then for any T 

ILdt + l|Vi.|r + \FEvf\,=T < CT(|T|lL,[„,rl.n) + l|G|lL(|„,r,xr))- 

Proof. Without loss of generality we can assume that F and G are both real. Multi¬ 
plying Dn = F with Vt and integrating over G we hnd 

^ (/ ^ j «( 7 U)^ + J G-fVt + j F Vt. (59) 

Using the Cauchy-Schwarz inequality on the second term on the RHS of fl5^ and 
recalling that a is strictly positive, we see that we can bound the hrst two terms by a 
multiple of /p G^. The other term on the RHS of (j5^ is bounded by -|- f^{vtY), 
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and the result then follows from Gronwall’s inequality (see, e.g., [241 §7-2.3]), using the 
fact that 6 > 0. □ 

In the proof of Theorem 11.81 below, the crucial microlocal ingredient will be the esti¬ 
mates on the wave equation with impedance boundary condition obtained by Bardos- 
Lebeau-Rauch [7]. These estimates involve a key geometric hypothesis, which is that 
every generalized bicharacteristic in the sense of Melrose-Sjostrand [55] eventually hits 
the boundary (or, in the more general setting of [7], the control region) at a point that 
is nondiffractive as dehned in [3, p.l037]B In our simple case of compact Euclidean 
domains, we remark that these hypotheses are always satished: 

Lemma 5.3. If G- C M” is a compact domain with smooth boundary, then every 
generalized bicharacteristic eventually hits the boundary at a nondiffractive point. 

Proof. We hrst observe that a generalized bicharacteristic in a compact Euclidean do¬ 
main must eventually change momentum. Adopting the notation of Hormander [361 
Dehnition 24.3.7], we claim that the only way the momentum can change along a gen¬ 
eralized bicharacteristic is when it hits the boundary at a point in "H U ^ \ Qd- Here PL 
denotes the “hyperbolic points” at which there is transverse reflection from the bound¬ 
ary, while Q\Qd denotes the set of glancing points that are not diffractive. To prove 
this assertion, we note that in the interior and at diffractive points (which together con¬ 
stitute the remaining parts of the characteristic set), we have 7 '(t) = Hpf^it)), where 
7 denotes the bicharacteristic and Hp the Hamilton vector held, which in this case is 
the constant vector held ^ • cla, in T*M"' (cf. Chapter 24 of [36]b 

Now we further note that on ^ \ Qd^ we have 7 '(t) = Hpipyit)) with Hp the “gliding 
vector held” of Dehnition 24.3.6 in [36] • This vector held still agrees with Hp unless 
’jit) G the points where contact with the boundary is exactly second-order. On the 
other hand, the “gliding points,” Qg = \ Qd, are nondihractive by the dehnition of 

Bardos-Lebeau-Rauch, since the second derivative of the boundary dehning function 
is strictly negative along the how at such points (cf. Dehnition 24.3.2 of [36]). □ 

Proof of Theorem \l.tA We begin by dealing with the case when a is positive. By [7] , if 
V satishes fl58|) with initial data in the energy space, then all energy norms of v enjoy 
exponential decay as t —)■ oo. Indeed, [T] Theorem 5.5 and Proposition 5.3] prove this 
result for the case when b is nonnegative, and then the result for 6 = 0 follows from [3 
Theorem 5.6], but we emphasize that in this latter case it is just the energy norm 

IIdII^ + ||Vnf + 

that converges to zero, while the value of the solution may converge to a nonzero 
constant, since this norm does not in general control the norm. 


^Note that the negation of “nondiffractive” in this sense is not the same as “diffractive” in the sense 


of 
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We let Tk denote the (unique) solution to the wave equation on R x satisfying 

(60a) 

[dn + 07(9* + h^)v^ = e“'^V(^)fi', (60b) 

nK(t,x) = 0, t < 0. (60c) 


where <yC is a cutoff compactly supported in (0,1) with f ip = 1. Then the standard 
energy estimate proved in Lemma 15.21 yields 

Il(t-4,|t + II V«.f + L., < ll/lli.,„, + llsIlLir, ■ 

Now since satishes the homogeneous wave equation for f > 1 with initial data at 
t = 1 controlled as above, [71 Theorem 5.5] yields, for some (5 > 0, 

ll(w)t||^ + llVn^ll^ + < Ce-^^( ||/||i2(f*) + ||^|li2(r)), t > 0. (61) 

Fourier transforming (l60l) gives 

(A + = -(p{k - K)f, 

{dn - ika'y + 67 ) = ^{k - K)g. 


Since 


U^v\ 


Li ^ 


the exponential decay estimate (1^ implies that 


||VJ-V|| + |fc|llV^“^.|| < ll/ll + lliill, k e R; 

here we have made no use of the boundary term on the LHS of fl6T|) . If the stronger 
Assumption 11.71 holds, we employ the more precise version of our Fourier transformed 
estimates: 

llvrt.ll' + Ill'll< Wff + \\gf, k e R. 

By the PoincarAWirtinger inequalitjH and the positivity of b, the left side controls 

||V^-V|r + (A:)'||rt.lf 

even at fc = 0, giving ns the stronger estimate (cf. discussion on pps. 1060-1061 of [7]): 

llV^-'n^ll + {k)\\j^-W\\ < ll/ll + lli^ll, keR. 

Taking k = k makes u = Vk the solution of the IIP (jHl) and yields the asserted estimate 
flTUD when a is strictly positive. This concludes the proof for a strictly positive. 

When a is strictly negative, the sign convention of the Fourier transform and the 
signs of the exponents in fl60D can both be changed to give the correspond estimate. 
(Alternatively, by taking the complex conjugate of the BVP ([9]), we can prove (fTOj) 
when the boundary condition 


{dn + iA:a 7 + b'y)u = g 


^Note that in employing the Poincare-Wirtinger inequality, we may estimate the average value of u 
by a multiple of ||Vu|| + |17m|| by writing it as a multiple of f uV ■ xdx and integrating by parts. 
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is imposed. If a is strictly negative, then we apply the bound above with a replaced by 
—a, and this yields the desired result.) □ 

We now prove Corollary 11.91 regarding the impedance-to-Dirichlet map, by using 
Theorem 11.81 in conjunction with a simple energy estimate. 


Proof of Corollary \1.9[ Reiterating the integration by parts used to obtain Lemma [5.11 
but now including the inhomogeneities, we hnd that applying Cauchy-Schwarz to our 
expression for the imaginary part of Jq fu yields for /c G M 


fc||\/a7M||^2(r) 

Now applying Cauchy-Schwarz and the estimates of Theorem 11.81 to the resulting ||m|| 
term on the RHS gives the desired estimate on /c^|| 7 m||^. 

The corresponding estimate for Vr(7M) follows from the analogous estimate to Lemma lT3l iii 
for bounded domains (the same proof employed by Necas applies). 

If b is strictly positive, we obtain the stronger estimate at A; = 0 by examining the 
real rather than the imaginary part of / u to estimate . □ 


91U 

+ 

/ fu 

Jr 


Jn 


Proof of Corollary \1.1(A We follow the argument in, e.g., [23l Theorem 2.5], [151 text 
between (3.3) and (3.4)]. The variational formulation of the IIP is 

hnd u G such that a{u,v) = F{v) for all v G iL^(fl), (62) 


where 


and 


a{u,v) := / Vu ■ Vv — k^uv — ik / ayuyn-l- / h'yu'yv, 


(63) 


F{v) ■- (/, v)n + {g, 7 v)r, (64) 

where (•, ■}q and (•, •)r denote the duality pairings on and T, respectively. Dehne the 
sesquillinear form ao(-, •) by 


ao(u,v) := / Vu ■ Vv + k^uv — ik / aynyn-l- / b'yu'fv. (65) 

Jn Jr Jv 

Furthermore, dehne uq G H^{VL) as the solution of the variational problem ao(now) = 
F{v) for all v G and dehne w G H^{VL) as the solution of the variational problem 

a{w,v) = 2k'^ J^uqv for all v G These dehnitions imply that the solution of 

(l6^ satishes u = uq + w. 

Since b is nonnegative, Reao(n,n) = I1 u||^i(q); thus, by the Lax-Milgram lemma, 
ll'*^o|li)fi(f 7 ) dehnition of w implies that w satishes the IIP with 

g = 0 and / = 2k‘^uo, and thus the bound fiTOj) implies that ||w||// 1 (q) ^ II'*^o|Il 2 (q). 

Combining these bounds on uq and w, we obtain 


W'^WhUu) ~ 1^1 ll-^ll(irRo))' • (66) 

The result on the inf-sup constant ([T3|) then follows from, e.g., [TOl Theorem 2.1.44]. 
The bound flT^ follows from fl66]) using the dehnition of F fl6T|) . □ 
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Proof of Corollary \1.11{ The bound ffldl) follows from combining the bounds flTOll and 

where the latter is proved in, e.g., [33l Theorem 2.3.3.2, page 106]. □ 

We now impose the homogeneous impedance boundary condition, and consider the 
operator Ri^j^{k) : L‘^{PL) —)■ L‘^{VL) dehned by Ri^ri{k)f = u where u is the solution to 
(A + k‘^)u = f satisfying — iri'j)u = 0. 

Following the discussion in 1|I1 we now proceed with the assumptions that a and b 
are both strictly positive (i.e. fl5Tl) and fl55ll hold), so that Ryriik) is well dehned for all 
Im /c > 0. 

We break the proof of Theorem 11.121 down into several steps; the hrst step is to prove 
that Rj^r]{k) is holomorphic on Imk > 0. 

Lemma 5.4 (Analyticity for Imk > 0). The operator family Ri^r^{k) : L‘^{Q) 
L‘^{Q) with boundary condition 

dnU — i{ka + ib)'yu = 0, (67) 

where a, b are real-valued C°° functions with a strictly positive on T and b nonnegative, 
is holomorphic on \mk > 0. 

Proof. First note that the standard variational formulation of the IIP satishes a Carding 
inequality. Indeed, the sesquilinear form is given by fl63|) and so, since b is non-negative 
and Im fc > 0, we have 

Rea(n,n) + (1 + k'^) 11^11^2(0) > Ikll^i(o) 

(note that we are using the unweighted norm on H^{Q) since we are allowing for k 
to be equal to zero). Fredholm theory then gives us well-posedness of the BVP as a 
consequence of the uniqueness result in Lemma EH] (see, e.g., [501 Theorems 2.33, 2.34]). 
Analyticity follows by applying the Canchy-Riemann operator d/dk to the equations 
(A -|- k‘^)u = / and dnU — \{ka + ib)'-^u = 0 : we hnd that du/dk must satisfy the IIP 
with zero interior and boundary data, hence by the uniqueness proved above, it must 
vanish. □ 

We now use a simple perturbation argument to get the existence of a pole-free strip 
beneath the real axis. 

Proof of Theorem \1.12[ Lemma 15.41 states that Ryn{k) is holomorphic on \mk > 0, 
while Theorem 11.81 yields the estimate (1T5|) for all fc G M (crucially using Assump¬ 
tion [FT]). We can now perturb using this estimate to extend to analyticity below the 
real axis, but we will need to consider the full inverse map on both interior and bound¬ 
ary data (and in so doing, will in fact prove a stronger result than stated, involving 
both interior and bonndary inhomogeneities). For the (unique) solution of the IIP 

(A -h k‘^)u = /, {dn - ifcay -F b^)u = g 



we set 


u 

'JU 
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Then Corollary 11.91 shows that for /c G M, 

Ri^r^ik) ; L^{Q) © L\r) ^ Hl{Q) © Hl(r). 

Now for z ^ C we may try to solve 

(A + (A; + zY)u = f, {dn - i{k + z)a'y + b'j)u = g 
by perturbation; we easily see that this is equivalent to 

(A + k'^)u = f — {2kz + z^)u, {dn — ika'j + b'j)u = g + iza^ju. 
Hence, applying Rj^riik), we wish to solve 


u 

7M 


= Rl,r^{k) 


f — {2kz + z‘^)u 


= Ri,n{k) 


Rj,,{k)M{z) 


u 

'yu 


g + iza'ju 

where 

'2kz + z^ 0 
0 —iza 

We can solve this by Neumann series (hence for a holomorphic solution with the same 
A;-dependent estimates as on the real axis) so long as, say. 


M(z) = 


Rj,^(k)M(z) 


-1 


< 1 / 2 . 


Since Rj^^{k) has norm bounded by C{k) ^ on © L^, this requires only that \z\ < e 
for some e > 0. Restricting to the case g = 0 gives the stated result. □ 


Lemma 5.5 (Sharpness of (ITOl) when / = 0 and hi is a ball). In for any d> 2 

there exist families of solutions u to the interior impedance problem in the unit ball 
with boundary inhomogeneity g : 

Am + k'^u = 0 in and dnU — vg'yu = g on (68) 

with 

^ll'^llL2(B<i) ~ llfi'llL2(Sd-l)- 

Proof. Fix any spherical harmonic (^(6^) on with eigenvalue —/i^. Then the function 

u{r,6) = Jy{kr)ip{6) 
solves the Helmholtz equation in B‘^ if we set 

z/ = ^V(d- 2)2 + 4/i2. 

We will let A: —)■ cxo while letting /i (and hence u) remain hxed. 

The function u thus satishes the HP (with rj = k) where 

g = {dr — iA:)M|r=i. 

Now we let A; —)> oo and examine the asymptotics of u and g. Since (see, e.g., m 
Equation 10.17.3] for the standard Bessel function asymptotics employed here) 


u = (p{9)r^ Y ^-( coscu + 0((rA:) ^)) 
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with 


we have 


U = rk — -UTT --TT 

2 4 


> (69) 

a.s k ^ oo with u hxed. On the other hand, using the asymptotic expansion of J(, yields 


|W|lL2 


hence at r = 1 we have 


drU =—(p{6)r^ sincj + 0(/c ^)), 

{dr — ik)u ~ ^{9)\j — ( coscjo + i sinwo) 


with uq = k — ^UTT — Thus, 


II - i/c)M||^2(5d-i) ~ 

Comparing to fl69|) yields the desired estimate. 


□ 


Remark 5.6 (Extension to inhomogeneous problems). The results of this sec¬ 
tion hold equally well, with identical proofs, if we generalize the fiat Laplacian to an 
inhomogeneous and/or anisotropic operator with smooth coefficients such as 


'^dia''^{x)dj, 

with a®'^ (x) strictly positive definite. The only difference is that we then need to impose 
an auxiliary geometric hypothesis, as Lemma 15.31 no longer applies. In this setting, 
motion along straight lines is replaced by the Hamiltonian dynamical system 


Xi{t) = 

iiit) = 


1 da^/x) 

9 / . fkfl- 


dxi 


(70) 


It may easily be the case that trajectories of this system—which are lifts to the cotangent 
bundle of geodesics with respect to the Riemannian metric a^/x )— fail to reach the 
boundary at a nondiffractive point or indeed at all (e.g. may be locally isometric 
in some region to more than half of a round sphere). Thus, we simply need to impose 
geometric control by the boundary as a hypothesis: we insist that all trajectories of (|7ni) 
do reach the boundary at a nondiffractive point. The rest of our results then follow as 
in the flat case. 


6. Boundary integral equations for the exterior Dirichlet and 

Neumann problems 

In this section we derive both the integral equation flTTD for the solution of the exterior 
Dirichlet problem and the analogous equation for the solution of the exterior Neumann 
problem. We then give a new proof of the decomposition (ITSD (which is more intuitive 
than the proof in [l3]), and we then prove an analogous decomposition for the integral 
equation for the Neumann problem. 
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We note that there are now many good texts discussing the theory of integral equa¬ 
tions for the Helmholtz equation, for example [50], [70], m, |3H]; we will use na as 
a default reference (since it, like us, is concerned with the high-frequency behaviour of 
these integral operators). 

If M is a solution of the homogeneous Helmholtz equation in 11+ then an application 
of Green’s formula yields 

u{x) = - J^^k{x,y)d:^u{y)ds{y) + J^^^^^j+u{y)ds{y), x e H+, (71) 

(see, e.g., [131 Theorem 2.21]), where ^k{.x,y) is the fundamental solution of the 
Helmholtz equation given by 

^k{x,y) = ^H^^\k\x-y\), d = 2, (x, y) = d = 3. (72) 

Taking the exterior Dirichlet and Neumann traces of (ITT]) on T and using the jump 
relations for the single- and double-layer potentials (see, e.g., [131 Equation 2.41] we 
obtain the following two integral equations 

7+w (73) 

and 

{^-I + D'^dtu = Hki+u, (74) 

where S'^, are the single- and double-layer operators, is the adjoint double¬ 
layer operator, and 77^ is the hypersingular operator. These four integral operators are 
defined for 0 e T^(r), ip G 77^(r), and a; G T by 

Ski^ix) := j d^k{x,y)ip{y)ds{y), DkPix) := j (75) 

I ^ I 

When T is Lipschitz, the integrals dehning Dk and D'f. must be understood as Cauchy 
principal value integrals and even when T is smooth there are subtleties in defining Hk'ip 
for Ip G T^(r) which we ignore here (see, e.g., [HI §2.3]). 

6.1. The Dirichlet problem. In the case of the Dirichlet problem, the integral equa¬ 
tions fl73]) and fl74p are both integral equations for the unknown Neumann trace dp^u. 
However fl73l) is not uniquely solvable when —/c^ is a Dirichlet eigenvalue of the Lapla- 
cian in D_, and fl74p is not uniquely solvable when —k'^ is a Neumann eigenvalue of the 
Laplacian in D_. (This is because if w solves the interior Helmholtz equation. Green’s 
formula yields 

Qj + = Skdpw; 

hence existence of nullspace of these operators is equivalent to existence of Dirich- 
let/Neumann eigenvalues.) 
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The standard way to resolve this difficulty is to take a linear combination of the two 
equations, which yields the integral equation 

= Bk,r,7+U 

( 77 ) 

where 


■=\l + Dk- ivSk 

( 78 ) 

and 


Bk,q ■ Hk T iy ^2^ ^k^ • 

( 79 ) 

If y G M \ {0} then the integral operator A'^ ^ is invertible 

(on appropriate Sobolev 


spaces) and so ffT7)) can then be used to solve the exterior Dirichlet problem for all 
(real) k. Furthermore one can then show that if 17 G M \ {0} then is a bounded 
invertible operator from H^{T) to itself for —1 < s < 0; [T^ Theorem 2.27]. 

For the general exterior Dirichlet problem it is natural to pose Dirichlet data in 
7/V2(r) 

(since 7 +M G i7^/^(F)). The mapping properties of and (see [121 The¬ 
orems 2.17, 2.18]) imply that : FF^+^(F) H^{r) for —1 < s < 0, and thus 

Bk^n7+u G i7“^/^(F). This indicates that we should consider fl77|) as an equation in 
77-i/2(r). 

Unfortunately evaluating the inner product numerically is expensive, and 

thus it is not practical to implement the Galerkin method on fll7p as an equation in 
7 ^-i/ 2 (f) ("fQj. short overview of proposed solutions to this problem, see [131 §2-11]) 
Fortunately, we can bypass this problem in the case of plane-wave or point-source 
scattering. Indeed, in this case y+u G H^{T) and G T^(F) [131 Theorem 2.12]. 
Since and A'^ .^d'^u are then in T^(F), we can consider fl77|) as an equation in 

L^(F), which is a natural space for implementing the Galerkin method. 


6.2. The Neumann problem. In the case of the Neumann problem we can view 
fl77|) as an equation to be solved for y+n. Indeed, given d^u G we have 

G 77“^/^(F) and Bk^^il+u G 77“^/^(F). The equation d77|) can then be cast as 
the variational problem on 77^/^ (F): hnd (j) G 77^/^ (F) such that 

{Bk,q(l), '0)r = ^|J)^ for all ^|J G 

where recall that (•, •)r is the duality pairing between 77“*(F) and 77*(F) for 0 < s < 1. 

Although this gives a practically-realizable Galerkin method, the fact that Bk,q is a 
hrst-kind operator means that the condition number of the discretized system depends 
on the discretization and thus it is desirable to precondition the equation with an 
operator of opposite order before discretizing (see, e.g., [76l §13] for a discussion of this 
technique in general). 

For Bk f^ this strategy amounts to multiplying fl72|) by an operator R ; 77“^(F) 
L^(F) and then adding it to —iy multiplied by fl73]) . This results in the equation 




( 80 ) 
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where 

Bk,ri := RHk + ir] Qj - 

and 

K, ■= R (h + fld - i-nS,. 

The mapping properties of R and the boundary integral operators Sk,Dk,D[, imply 
that both Bk^r) and are bounded operators mapping T^(r) to itself, and thus, in 
the case when d^u G T^(r), flSU]) can be considered as an integral equation in T^(r). 
Of course, R must satisfy some additional conditions to ensure that flHOl) has a unique 
solution for all k > 0. 

The most common choice is to take R = So, motivated by the Calderon identity 

SoHo = ~I + Dl 

f [T3l Equation 2.56]) and the fact that So{Hk — Hq) is compact (since — Hq has a 
weakly singular kernel; see m Equation 2.25]). 

The choice R = Sik was proposed in jin], and further used and analyzed in, e.g., 
0. 0. Other choices for R include principal symbols of certain pseudodifferential 
operators [9|, and (for the indirect analogue of flHOji l approximations of the NtD map 
0§8|. 

6.3. Decompositions of inverses of combined potential operators. The decom¬ 
position (lT8l) of in terms of P^tN implicit in much of the work 

on but was hrst written down explicitly in [HI Theorem 2.33], along with the 
analogous decomposition for B^]^ (as a special case of the decomposition of the inverse 
of the integral operator for the exterior impedance problem). 

In Lemma [6. II below we provide an alternative, more intuitive, proof of these decom¬ 
positions. We also give the analogous decomposition of the operator B^]^ in terms of 

and where the operator ■ L^(r) —)■ T^(r) maps g G T^(r) to the 

Dirichlet trace of the solution of the BVP 

Au + k‘^u = 0 in r2_, Rd~u — ig'y^u = g on T 

(assuming appropriate conditions on R are imposed so that this BVP has a unique 
solution for all k > 0). 

Lemma 6.1. ITe have the following expressions for the inverses of combined-potential 
operators: 

i^k,ri) ^ “ i^DiN ~ ''^'h)PltDi 

(St,,)-' = Pi^o -a- ilPmD)PnD, (82) 

(ft,,)-' = PS,dR-' -(I- i’lPS,DR-')Pr:D'‘- (83) 


Proof of Lemma \6.1[ We recall (e.g. from Section 2.5 of [13]) the formula for the interior 
and exterior Calderon projectors, which project onto pairs of Dirichlet and Neumann 
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data for solutions to the Helmholtz equation in and (with radiation condition) 
respectively. In terms of layer potentials, we may write these operators as 

n, = i/±A4. 

(Here we have departed from the notation of [13] for the Calderon projectors—these 
authors use P ±—as the letter P is somewhat overloaded.) 

These dehnitions imply that 

{-ir] 1) ■ n_ = J . 

Hence 

(-ir, 1) ■ n_ Cj ^ = g. (84) 

On the other hand, since n„ projects to Cauchy data for the interior Helmholtz 
problem, we assuredly hnd that 




means that 



are Cauchy data for the interior impedance problem, hence we may rewrite 


n_ 



(PuSi9)\ 

W’MJ ■ 


Since n+ -I- n_ = /, we now hnd that 


(85) 


n 


+ 



(4>-PnSi9)\ 

W - PiMJ • 


Note that the RHS is now guaranteed to be Cauchy data for a solution of the exte¬ 
rior Helmholtz equation (with radiation condition) and hence we may write its two 
components in terms of one another via the maps and P^^d- 

Now we split into the special cases of 0 = 0 or = 0. In the former case we have 


H 


+ 


-PnSi9) 


'PutN^P 110(9)) 


(where we have written the second component in terms of the hrst using Pdin)- Thus 


= [-IT] 

= (-ih 
= (-ih 


■)■(:) 

i).(n, + n_)(°) 


!)• 


-pp:]{9) 


-P: 


ItD'' 
(p-,V 
DtN\^ItD 


(9)) 


+ 9 
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where we have used 


Thus 


to evaluate the n__ term. Likewise, when '0 = 0 we have 

\9) 


n. 


PuSi 


DtN 


(0 


^ItD 


(9)) 


= {-if] 1) • 

= {-iv i)-(n+ + n„)(^^^ 


= 1) • ( 


p+ 

^DtN 


^ PltDi9) 

- PnSi9\ 

In both cases, solving for -0 (respectively 0) and recalling fIMl) gives the desired expres¬ 
sion in terms of g (in the latter case, we use that 0 = o P^^^0). 


Finally, to obtain the formula for we apply the same argument as for Bf^ but 


where we consider 
throughout, rather than 


(—ir; i?) ■ n_ 
—irj l) ■ n„. 


□ 


The estimate ^ analogous to the estimate flT^ on (A^ ^ is as follows. 


Lemma 6.2. Let r2_|_ C d 
that 7] satisfies Assum.ption \l.fA 


= 2, 3, be a smooth, nontrapping domain and suppose 
Then, given ko > 0, 


B 


-11 




< 


\k\ 


1-/3 


( 86 ) 


for all \k\ > kfj, where fi is as in Theorem \l.,5[ 


Since this integral operator is not used in practice, however (as explained in §6.2p . 
we do not include the proof. Note that an estimate from to H^P{T) can be 

obtained from (1861) by interpolation. 

The decomposition of given by fl83l) below and the sharp bounds on in 

Theorem 11.51 reduce the problem of bounding ||50,^||L2(r)^.L2(r) to that of bounding 
Phd^ for the different choices of R, however we do not pursue this further here. 


7. Concluding remarks: the conditioning of 

In l ll.di we stated that the present paper combined with the recent work of Galkowski- 
Smith and Galkowski almost completes the study of the high frequency behaviour of 
11^4). ,^11 and ||(^p^)~^|l, and thus of the condition number 

cond(A'^^^) := ||^fc,r,IUFr)^LFr) IK^^pr?) ^\\ L '^ ir )-^ L ^{ r )- (87) 

We conclude this paper by justifying this remark in 1 17.11 but then also questioning in 
1 17.2l whether the condition number is an appropriate object to study in relation to 
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7.1. Upper bounds on cond(74'^^). We begin by recalling the recent sharp bounds 
on ||*S'fc||i:, 2 (r)^ 2 , 2 (r) and \\Dk\\L'^{T)^L'^{T) proved in [2H1 Theorem 2], [HU Theorem A.l]. 
(Note that \\Dk\\L‘^{T)^L‘^(T) = ||-D^||L 2 (r)^L 2 (r), and so these bounds are sufficient to 
bound ||Wfc jL 2 (r)^L 2 (r).) 

Theorem 7.1. ([28l Theorem 1.2], [HU Theorem A. 1], m Theorem 4.4]) With 
and T defined in §7.71 ifT is a finite union of compact embedded hypersurfaces then 
there exists ko such that, for k > ko, 

ll‘^fcllL 2 (r)-j.L 2 (p) ^ A; ^ log/c, Il-Dfc||p 2 ('p)_^p 2 ('p) ^ logfc. 

If V is a finite union of compact subsets of C°° hypersurfaces with strictly positive 
curvature, then 

ll‘^fcllL 2 (p)_>.p 2 (p) A; ! logfc, ||-Dfc|li, 2 (p)_^p 2 (p) ^k ^ \ogk. 

Moreover, modulo the factor log k, all of the estimates are sharp. 

Note that in 2-d the sharp bound ||S'fc|]p 2 (p)_j.p 2 (p) < proved in [121 Theorem 

3.3]. 

Combining these bounds with the bounds on || (|2T|) and ffT9|) . as well as 
bounds when T is the circle or sphere obtained by ra. iii. 0 (see the review in [131 
§5.4]) we obtain the following theorem. 


Theorem 7.2 (Upper bounds on the condition number). 

(a) Let be star-shaped with respect to a ball, with T piecewise smooth. When d = 2, 
*/ 

/c^/^ogfc <\r]\<k 


then 


When d = 3, if 


cond(A'fc,,) < /c^/2. 
^\v\^k 


( 88 ) 


then 

cond(Afc ,^) < log k. (89) 

(b) IfLl- is nontrapping and rj satisfies Assum.ption \1.6\ (which includes the case \ri\ ~ 
k), then fl88|) holds when d = 2 and holds when d = 3. 

(c) If fl- is star-shaped with respect to a ball, T is the finite union of smooth surfaces 
with strictly positive curvature, and 


< LI < 


then 

cond(Afc ,j) < log k. (90) 

In particular, iffl- is a 2- or 3-d ball (i.e., T is the circle or sphere) then cond(A'^_^) < 
k^^^ when 


k^'^ < \v\ ^ k. 
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Earlier we stated that this theorem “almost completes” the study of cond(A^ ^). One 
thing that is missing is a lower bound on cond(y4'^,^) that shows the choice |? 7 | ~ fc is 
optimal. Indeed, in 2-d, if E contains a straight line segment, then by [T^ Theorem 
4.2] 

lU' II + 

II fc.^?llL2(r)^.L2(r) ^1/2 ^ I ^ / 

as k ^ oo, uniformly in \ri\. The only existing lower bound on ||(^fe^^)~^|| is ||(^fc,^)~^|| > 
2, which holds if a part of T is [T2l Lemma 4.1], and with this alone we cannot rule 
out the possibility that cond(y4';.^^) for a choice of Ir^l k but > k^/'^logk 

(although we do not expect this to be the case). 

7.2. Should we really be interested in the condition number? To be concrete, 
we consider solving numerically the integral equation fll7p (as an equation in L^(r)) 
via the Galerkin method, i.e. given a sequence of hnite-dimensional nested subspaces 
V/v C T^(r), we seek E Vn such that 

{A'f^ r,VN-,WN)L^{V) = {fk,'q-,W n)L'^(T) for all Wn G Vn- (91) 

We restrict attention to the case when Vat consists of piecewise polynomials (and so 
we do not consider, e.g., subspaces involving oscillatory basis functions; see, e.g., [13] 
and the references therein), and furthermore we only consider the h-boundary element 
method (BEM) (i.e., the piecewise polynomials have hxed degree but decreasing mesh 
width h). 

Given a basis of V/v, equation (1^ becomes a system of linear equations; for simplicity 
we do not consider preconditioning this system. 

For the high-frequency numerical analysis of this situation, there are now, roughly 
speaking, two tasks: 

(1) We expect that the subspace dimension N (~ must grow with k in 

order to maintain accuracy, and we would like k- and ? 7 -explicit bounds on the 
required growth. 

(2) One usually solves the linear system with an iterative solver such as the gener¬ 
alized minimal residual method (GMRES); we expect the number of iterations 
required to achieve a prescribed accuracy to increase with /c, and we would like 
k- and r^-explicit bounds on this growth. 

Regarding 1: The analysis in [31] shows that there exists a C > 0 such that if 

h (^ll-Dfc||i2(r)-^//i(r) + 1^1 ll-^fcll ) ^IU2(r)^L2(r) < C 

then the sequence of Galerkin solutions is quasioptimal (with the constant of qua- 
sioptimality independent of k), i.e., 

pin - || 8 > - ; 

see [311 Corollary 4.1]. Therefore, minimizing 

(ll-^fcllL2(r)^.Ri(r) + 1^1 ll>^fell L2(r)^rji(r) ) \\i^'k,ri) ^IU2(r)^L2(r) (92) 

gives the least restrictive condition on h. 
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This is not quite the same as minimizing the condition number, but if we believe 
that the —)■ iJ^-norms of and Sk are proportional to the —)■ L^-norms (with 
the same constant of proportionality), as they are in the case of the circle and sphere 
at least (with “constant” of proportionality k), then minimizing fl9^ is equivalent to 
minimizing the condition number^ 

Two remarks: 

• In [31] bounds on the i/^-norms are obtained and it is shown that, if \ri\ ~ 

k and hi- is both and star-shaped with respect to a ball, then the quantity 
in (l9^ is bounded by in 2-d, yielding the condition for quasioptimality 
^^3/2 < gf ^pg circle/sphere, better bounds on the norms can be 

used to obtain the condition for quasioptimality hk^^^ < 1. In practice, one sees 
that the h-BEM is quasi-optimal when hfc < 1 (i.e., it does not suffer from the 
pollution effect), see, e.g., [211 §5], but this observation has yet to be established 
rigorously. 

• Here we have only talked about the h-BEM; the hp-BEM (where the polynomial 
degree, p, is variable) is less sensitive to the value of rj and the norms of 
and (^fc,j)~^; see [16], [52] for more details. 

Regarding 2: In the discussion above we noted that, in practice, h/c < 1 is sufficient 
to ensure fc-independent quasioptimality of the Galerkin method. Since N ~ ^ 

this condition implies that, as k increases, the size of the linear system must grow like 
k^-^) to maintain accuracy. Iterative methods, such as GMRES, are then the methods 
of choice for solving such large linear systems. 

For Hermitian matrices there are well-known bounds on the number of iterations of 
the conjugate gradient method in terms of the condition number of the matrix [321 
Ghapter 3], and for normal matrices there are well-known bounds on the number of 
GMRES iterations in terms of the location of the eigenvalues (which can be rewritten 
in terms of the condition number) [671 Theorem 5], [661 Gorollary 6.33] (how satisfac¬ 
tory these bounds are is another question, but they exist). In contrast, for non-normal 
matrices it is not at all clear that the condition number tells us anything about the be¬ 
haviour of GMRES (at least, there do not exist any bounds on the number of iterations 
in terms of the condition number of non-normal matrices). 

As a partial illustration of this in the context of Helmholtz integral equations, the 
recent work of Marburg BH. liai has emphasized that, at least for certain collocation 
discretizations of the integral equation fl77|) . used as an integral equation for the Neu¬ 
mann problem, the sign of rj affects the number of GMRES iterations (with rj = k 
leading to much smaller iteration counts that rj = —k). An analogous effect occurs for 
similar collocation discretizations of the integral equation (177|) used as an equation to 
solve the Dirichlet problem, with the choice of rj = k much better than rj = —k [48] . 
In contrast, the condition number estimates in Theorem 17.21 are independent of the 
sign of ? 7 , suggesting that the condition number is not the right tool to investigate the 
behaviour of GMRES. 

®The methods used to prove the bounds in Theorem 17.11 also appear to be able to prove the corre¬ 
sponding — >■ bounds with an extra factor of k on the right-hand sides [55]; thus the proportion¬ 
ality discussed above would hold. 
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A concept that does give bounds on the number of GMRES iterations for non-normal 
matrices is coercivity. On the operator level (for A'f,^ on L^(r)), coercivity is the 
statement that there exists an ak,n > 0 such that 

l(^M‘^>>i2(r)l > I10lli2(r) all 0 G L^{r), 

and the matrix of the Galerkin method then inherits an analogous property (see, 
e.g., m Equation 1.20]). If is coercive, then the so-called Elman estimate for 
GMRES [21], [20] Theorem 3.3], [Ml §1-3.2] can be used to prove a bound on the number 
of GMRES iterations required to achieve a prescribed accuracy, with the bound given 
in terms of and || ,^||^ 2 (r)^L 2 (r); see [H Equation 1.21]. 

It is not clear whether bounds on the number of GMRES iterations obtained via this 
method are sharp, and so far ^ has only been proved to be coercive when rj > k and 
is strictly convex (and under additional smoothness assumptions on T), so we do not 
yet know enough to make a provably-optimal choice of rj via this approach. However, 
we do know that the sign of rj does matter for coercivity. Indeed, when is a ball, 
A'j^.^ is coercive when rj = k [TS], but not when 7] = —k [TH §1.2]. The dependence 
of coercivity on the sign of r] is consistent, therefore, with the results of Marburg that 
indicate that the number of GMRES iterations for A'^ ^ depends on the sign of rj. 
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